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Abstract 

Experimental correlation between the rms nucleus charge radius and the neutron-proton asym- 
metry is discussed. Simple attempt of explanation of this correlation using a semi-empirical liquid 
drop model which takes into account nuclear compressibility, nuclear deformation and neutron skin 
effects is presented. It seems that such procedure could be used as a tool for determination of the 
nuclear equation of state. 
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The nuclear equation of state (EOS) and the nuclear reaction kinetics determine stellar 
structure and evolution. In particular, the EOS allows to estimate the critical parameters 
at which the collapse process is halted and the shock wave in a core collapse supernova 
becomes formed |2j]. But the equation of state is not only relevant in astrophysical 
aspects. Knowledge of EOS is a way to better understand the nuclear heavy-ion reactions 
0, 0] and some fundamental properties of heavy nuclei, like for instance binding energy or 
stability of neutron-rich nuclei p|. 

The equation of state for infinite, asymmetric nuclear matter usually represents the rela- 
tion between the energy density per nucleon (baryon) e and densities reached by neutrons 
and protons in the nuclear matter. In this relation the Coulomb energy is not taken into 

Q 

account. Many theoretical approaches [6] have shown that the properties of infinite asym- 
metric nuclear matter, at relatively low temperature, can be described by an approximate 
form of EOS: 

e(p, 5) = e(p, 5 = 0) + 5 2 e sym (p) (1) 

where the baryon density p is a sum of neutron p n and proton p p density. Isospin asymmetry 
5, is defined as 

g Pn Pp Pn Pp /g\ 

Pn + Pp P 

First term in the eq. (QJ) represents energy density per baryon associated to symmetric 
nuclear matter (5 = 0). This energy can be expanded around normal nuclear matter density 
Po as 

e(p,0)=e( Po ,0) + ^(^)V. (3) 

In the formula above Ks=o is the nuclear compressibility for symmetric matter, defined as 

2 d 2 e{p, 0) 

A 5=0 = 9p dp2 Ip=po • (4) 

Second term in eq. (pj) stands for the symmetry energy and can be also expanded around 
Po as 

( \ ( \ i L (P~Po\ , K sym (p-PoY . 

e sym (p) = e sym (Po) + 3 I — J + ( — J + - (5) 
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where variables L (slope) and K syrn (curvature) characterizing the density dependence of 
the nuclear symmetry energy are given by 

T _ o de sym (p) 
~ — dp — lp=p0 ^ ' 

and 

t^sym — VPs q^ 2 \p=po \' ) 

The coefficients po, e(po,0), e sym (po), Ks=o, L, K sym allow us to determine the EOS 
around the normal nuclear density. Unfortunately, values of the slope L and the curvature 
K sym are neither well determined experimentally nor accurately estimated theoretically |f|. 
Data from the giant monopole resonances suggest that K sym varies from -566 ±1350 MeV 
to 34 ±159 MeV (jf. So far, experimental predictions for the L value are rather scarce. 
In the recent paper [a| the isospin transport calculations compared with the MSU diffusion 
data suggest L « 60. In theoretical models L varies from -50 up to 200 MeV 0|. Moreover, 

heoretical estimates for K sym are model-dependent and vary from about -700 to ±466 MeV 

lj 

In the present paper we will show a way to put some constraints for possible values of 
these parameters. In order to do it, we compare the experimentally measured correlation 
between root mean square radii of nuclei and isospin asymmetry, to the similar correlation 
resulting from a modified liquid-drop model (LDM). 

One of the most elementary properties of a nucleus is the spatial nucleon density distribu- 
tion. Unfortunately, this distribution cannot be obtained experimentally in a straightforward 
way. In the case of protons, knowledge of the Coulomb interaction allows to measure with a 
high precision the momenta of the charge distribution. From experimental data one usually 

1 /2 

fixes the root-mean-square (rms) radius (r^J . In this paper we will use a part of the set 
of 799 ground state nuclear charge radii presented in (ll|. This experimental data were 
obtained from the elastic electron scattering, muonic atom X-rays, K a isotope shift, and 
optical isotope shift. Errors associated with the data are usually much lower than 1%. For 
neutrons, the differences between neutron and proton density distributions at large nuclear 
radii in stable nuclei may be determined e.g. from antyprotonic atoms [lj| 

In the first column of Fig. [2 we plot a correlation between the experimentaly measured 
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(r^j ^ I \/~A and isospin asymmetry I = (N — Z) / (N + Z) , where N and Z denote neutron 
and proton numbers, respectively. In Fig. we plot these correlations for all available 
nuclei with masses A = N + Z > 100. We can notice here some branches which are roughly 
mutually parallel. These branches are much better visible on the next drawings in this 
column (Fig. [it, e, g) showing the correlations for nuclei with selected mass intervals. 

Let us try to explain the nature of the observed correlations. We consider here the 
influence of three effects: nuclear deformation, dependence of nucleon density on isospin, 
and neutron skin. 

At the beginning we will study the influence of nucleus deformation on the proton rms 
radius. In our analysis we neglect the diffuseness of the surface and we assume only ground- 
state Quadrupole deformation described by "/3 parametrization" (a spherical harmonic expan- 
sion) |l2|]. Such an aproach should be sufficient for rather small ground state deformations. 
Now we can write the proton mean-square radius in the form 



M 1 + (8) 



arameters j3(A,I) for considered set of nuclei are taken 
1^ | . Obtained correlations between deformation (3 2 and 



where Rq is the average radius. 3 
from model calculations used in 
isospin / are presented in the right column in Fig. QJ As we can see from the upper graph 
(Fig. |T]b) for almost all masses (3 2 is less than 0.1 . Formula (jHJ) suggests that for a given R , 
nuclear deformation does not modify the (r£) value more than 1% , however (see Figs 
IUl,f,h) it can enhance or suppress the observed correlations depending on the selected mass 
range. Even more importantly, it can also change the correlation sign (see e.g. FigHk and 

FigET)). 

One can suspect, that such a strong correlation observed in the first column could be 
associated with the form of the nuclear equation of state. To examine this possibility we 
assume that the considered correlation can be a result of differences in the density of nuclear 
matter in the central region of nuclei due to different values of the surface tension and 
Coulomb interaction. 

Saturation of the nuclear matter density in the central region of nuclei provides a base 
for description of the nuclear binding energy in the liquid-drop model. This model was 
formulated by Weizsacker (1935), Bethe (1936), and Bohr (1936). Up to now various mod- 
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ifications were introduced in this model and some of them consider also compressibility of 
nuclear matter fl^ . 

The liquid-drop model in its simplest version describes the nuclear binding energy B as 
a sum of the volume, surface and Coulomb terms. 

B = B V + B S + B c , (9) 
where B v is the dominant volume ingredient usually written as, 

B v = a v (l - k v I 2 )A, (10) 

and a v and k v are coefficients. The term k v I 2 takes into account the lowering of the binding 
energy caused by different number of protons and neutrons in the nucleus. 

Next term is (9) takes into account lower binding energy for nucleons located on the 
nuclear surface, and can be expressed as 

B s = a s (l - k s I 2 )A 2 ^f = [a s (l - k s I 2 )A^ 3 ] Af, (11) 

where a s and k s are coefficients. Factor / = (1 + |/? 2 — ^P 3 ) takes into account the nuclear 
deformation (see [3|). We can notice that the surface correction per one nucleon (square 
bracket in (fTTJ)) vanishes with the increasing A. 

Coulomb interaction between protons appears in the next correction, which is propor- 
tional to Z 2 and inversely proportional to the radius of the nucleus. This term can be written 
as 



Z 2 

gA (12) 



C A±I\ 

where a c is a parameter, and g — (1 — |/3 2 — y^/3 3 ) similarly like / in fTT]) describes possible 
influence of the nuclear deformation. Comparison of (jH)) to the experimental data (e.g. [l^| 
allows to determine coefficients presented in formulas (fTfl|) - (fT2|) as: 

a v = -15.677 MeV, a s = 18.560 MeV, a c = 0.717 MeV 

k v = k s = 1.79 MeV. (13) 
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It should be stressed, that LDM is dedicated to a finite drop of nuclear matter with a 
constant density. From the analysis of experimental data we can deduce that apart from 
nuclear surface the nucleus density is roughly constant and for mass A > 40 it has a Fermi- 
like shape distribution. The central region density itself somewhat varies for different nuclei. 

The form of the volume part of the binding energy fTUj) is in fact related to the EOS at 
the ground state density. For a v = e(po, 0), a v k v = e svm (po), p = po, and 5 2 = I 2 the energy 
obtained from EOS JT} is equivalent to the volume energy per nucleon in the LDM (fT0|) . For 
our consideration we reformulate the liquid-drop model in a way that allows to determine 
the central region density for a given nucleus. We want to obtain the central region density 
as a function of coefficients po, e (po>0), e sym (po), Kg =0 , L, K sym , and 5 = 1, determining 
the nuclear EOS. In principle neutron and proton density distributions may have different 
radii, R n , R p , what gives the so called skin effect. Lets assume for the beginning R n = R p . 

Modifying the LDM we assume that density in the central region is a function of global 
values N and Z (i.e. there is no local dependence). According to the idea of the LDM, we 
can say that the central region density p results mainly from nuclear interaction present in 
infinite nuclear matter, further being modified by the surface tension and Coulomb repulsion. 
Therefore, reformulating the liquid-drop model, we write the binding energy per nucleon 
b(p, I, A) in a form similar to 

b(p, I, A) = b v (p, I) + b s (p, I, A) + 6 c (p, /, A). (14) 
In the above formula the volume part of the binding energy is defined by EOS {TJ. 

b v (p,I)=e(p,5). (15) 

Let b v0 = e(po,5). Now, the correction for the surface effect we can rewrite in a new, 
density dependent form 

b s (p,I,A) = <(l-Kl 2 ) b ^(^) 2/3 A-^f 

b v o \p J 

I / \ 2 /3 

= a -rbMi)( P A A-Wf. (16) 

a v \ P J 

Factor [— \ ^ in (|TK|) is responsible for the surface change due to the change in the core 
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density. The corresponding coefficients are marked by prime to distinguish from parameters 
of the standard LDM. In principle their values can slightly differ from the original ones. 

The part describing the Coulomb interaction is modified only by factor ( according 
to the radius change induced by core density changes (Z = Z(I,A)). 

UpJ,A) = a' c ^(^\ ' 3- (17) 

Here the binding energy is determined by the density p and a new set of coefficients a' v , 
a' s , a' c , po, e(po, 0), e sym (po), K$=o, L, K sym . The density p for a given nucleus is found from a 
condition db(p,I,A)/dp \i,A=const= 0. In this way p and binding energy values are uniquely 
determined by a selected set of coefficients. Here, as a first approximation, we are taking the 
value of coefficients with the sign prime equal to corresponding coefficients without prime 
which are given by (TTTfl) . Coefficients po, e(po,0), e sym (po), are rather well known. In Fig. El 
we present the binding energy b(p,I,A), as a function of density, for nuclei with A = 150 
and three different Z values equal to 50, 60 and 75 (7 = 1/3, 7 = 1/5 and 1 = 0). For 
these test calculations we take Kg=o = 240 MeV and we check energy density behavior for 
all possible combinations of L= ±100 MeV and K sym = ±500 MeV what determine the 
symmetry energy (jUJ) . 

The saturation of nuclear density in the central region is characteristic for heavier nuclei. 
For light nuclei we can say that almost all nucleons are located on their surface. In order 
to omit the problem of high contribution coming from the surface we shall test formula (|T1*|) 
only for heavier nuclei, A>100 . 

Assuming _R ~ \J A/ p and using (jHJ we can rewrite relation between r rms and p in the 
form 




where a is a normalization constant. Formula (jTHjl neglects a neutron skin, which because 
of different radii of neutron and proton density distributions, surrounds the larger A nuclei. 
Lets define thickness of the neutron skin as 

d = R n - R p . (19) 

where R p = (r^ ^ , and R n = ^/f (O^is the effective proton and neutron radius, 
respectively. The central region nucleon density is now 
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N Z 

P = Pn+P P =|^ | + P ? (20) 

and instead of / = (N — Z)/(N + Z) one gets for the central region isospin I 

j_ N-Z(l+d/R p )\ ( 3Zd \ 

N + Z{l + d/R p f- y r (f) 1/3 A^) 

For numerical calculations we use r = 1.2049 (see 

As d seems to be much smaller than R p , one gets from fT9|) and ([20l 

/ 



R p 
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3A 

4:7ip 



Nd 



The Coulomb energy is inversely proportional to the radius and consequently b c in 1)1 7|l has 
to be written as: 



(22) 



Z 2 

b c (p, I, A) = a' c 





V3 / 


(?: 




V 



As an alternative for eq. (|T%J) one should use now: 



Nd . 

!+ ; ,i/3 \9- (23) 



r (e) V3 ^ 4 / 3 



4 n 

-7= = QJ J 



' Nd . 

V r (a) A 4 / 3 



Lets try to discuss the origin and slope of the observed correlation's as a function of 
different factors. 

In Figs. 01 and |U in columns marked "model 1", "model 2", and "model 3", we plot cor- 
relation's given by formula (J23J, for different sets of the EOS parameters. For calculations 
named "model 1" we neglect the neutron skin and deformations of nuclei (/3 = d = in 
formula (J2"4*j) ) whereas in "model 2" deformations are taken into account (/3 ^ 0, d — 0). 
In "model 3" we are trying to estimate the influence of the neutron skin (ft ^ 0, d ^ 0). 
Calculations were performed for the soft (-fC<y =0 = 120 MeV, red points) and for the hard 
(K 5=0 = 500 MeV, blue points) EOS. In Fig. E]L = -150 MeV and in Fig. HL = 60 MeV. 
As for nuclei p ~ po and in formula (5) the third term is a rather small correction, in both 
cases we take K sym = 500 MeV. 

1 /2 

For "model 1" the (rj) j\[A versus / correlation already contains a number of separated 
chains because of discrete values of the A and Z numbers. Nuclear deformations modify 
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slightly slopes of the observed chains and may result in additional splitting (see "model 2", 
180 < A < 210). As seen from the "model 3" calculations taking into account the neutron 
skin introduces quite important changes. 
Conclusions 

An experimental correlation between the nucleus rms charge radius and the isospin of 
heavier nuclei has been observed. It seems to be related to interactions characteristic for 
the nuclear matter equation of state modified by confinement inside of an atomic nucleus. 
These interactions are responsible for nuclear deformations and creation of the neutron skin. 
We propose explanation of the observed correlation on the basis of a modified nuclear liquid 
drop model. Inspection of studied correlation's (Figs El and 0J) suggests a possibility of 
evaluating the EOS parameters. A precise determination of these values requires a more 
precise search and a more complex analysis which will be subject of a future work. Such 
investigation probably requires an extension of the LDM by additional terms taking into 
account such corrections as for example pairing, shell effects and diffuseness of the nucleus 
surface. Although b s in formula (fTK|) decreases for large A values, it could be interesting to 
introduce here an additional correction related to the local isospin of the surface. 

Acknowledgments 

The author is indebted to P. Pawlowski for solving some computing problems and to J. 
Blocki, K. Grotowski, M. Misiaszek, and P. Pawlowski for careful reading of manuscript and 
helpful discussions. This work was supported by the Polish-French collaboration IN2P3, 
grant no. CSI8. 

Figure 1: Correlation between the neutron-proton asymmetry and (first column) the experimentally 
measured nuclear rms radius, or (second column) the calculated parameter of deformation. 

Figure 2: Binding energy per nucleon vs density p for different equations of state (see text). 

Figure 3: Comparison between the measured and calculated nuclear rms radius vs /, for L = 
-150MeV and K sym = 500MeV (see text for details). 
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Figure 4: Comparison between the measured and calculated nuclear rms radius vs /, for L = 60MeV 
and K sym = 500MeV (see text for details). 
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